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Spin precession has been used to measure the transmission time r over a distance L in a graphene 
sheet. Since conduction electrons in graphene have an energy-independent velocity v, one would 
expect t > L/v. Here we calculate that r < L/v at the Dirac point (= charge neutrality point) 
in a clean graphene sheet, and we interpret this result as a manifestation of the Hartman effect 
(apparent superluminality) known from optics. 
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I. INTRODUCTION 

The precession of the electron spin in a magnetic field 
provides a clock for the study of the electron dynamics 
[l[ . This so-called Larmor clock 0, 0| is a particularly 
useful tool in a quasi-two-dimensional system, when one 
can use a parallel magnetic field to avoid perturbing the 
dynamics by the Lorentz force. The single-atomic layer 
of carbon atoms known as graphene is the ultimate two- 
dimensional system Spin precession was used suc- 
cessfully by Van Wees and collaborators to measure the 
diffusion time through a disordered graphene sheet [H, Q . 
For a mean free path I small compared to the separation 
L of the source and detector contacts, the diffusion time 
to ~ tqL/1 is larger by a factor L/l than the ballistic 
time of flight tq ~ L/v, with v — 10 6 m/s the energy 
independent Fermi velocity in graphene. 

In a clean graphene sheet, when I 3> L, the diffusive 
dynamics becomes ballistic, at least for Fermi energies 
Ef away from the Dirac point (ef = 0). At the Dirac 
point the dynamics in graphene is called "pseudodiffu- 
sive" : conductivity and shot noise suggest diffusive trans- 
port even in the absence of any disorder Q . In this paper 
we theoretically address the question of what spin preces- 




FIG. 1: Schematic top view of a graphene sheet with four fer- 
romagnetic contacts numbered 1,2,3,4; arrows indicate the 
direction of magnetization. The ratio /31//41 of currents from 
contact 1 into contacts 3 and 4 measures the spin precession 
time t' 2 ' in an in-plane magnetic held B. An alternative ge- 
ometry, with the magnetization in contacts 3, 4 aligned per- 
pendicularly to the magnetization in contacts 1,2 (and still 
perpendicularly to B), measures the time t' 1 ^ through the 
ratio (hi - L a i)/(l4,i + I'm)- 



sion can tell us about the dynamics at the Dirac point. 

While the notion of pseudodiffusive dynamics might 
suggest a scaling r oc L 2 for the transmission time r 
at the Dirac point, such quadratic scaling is forbidden 
by dimensional arguments. In the absence of disorder 
there is only a single length scale L at ep = 0, sor = 
constant x L/v is the only quantity with dimensions of 
time. As we will show, the proportionality constant is 
< 1, so t < L/v — as if electrons could propagate at 
speeds > v. 

The optical analogue of this anomalously short trans- 
mission time, with v replaced by the speed of light, is 
called superluminality or the Hartman effect [1, 0] . As 
explained by Winful there is no violation of rela- 
tivity because the transmitted waves are not propagat- 
ing but evanescent. Graphene would offer an interesting 
possibility to observe this paradoxical effect in the solid 
state. 

In the next sections we formulate the scatteringprob- 
lem in a clean graphene sheet at the Dirac point [7], [ll| , 
and calculate the transmission time r measured in a 
weak-field spin precession experiment over a distance L. 
We then perform a separate calculation of the mode- 
dependent Wigner-Smith delay time r n , which is directly 
defined in terms of the scattering matrix [l2l [13 (with- 
out reference to spin precession). This is the quantity 
studied in the optical context. 

We demonstrate that r is the weighted average of r„, 
weighted with the mode-dependent transmission propa- 
bility T n . More precisely, depending on the relative 
alignment of the magnetization at the two ends of the 
graphene sheet, the precession experiment measures ei- 
ther r^ 1 ) or t( 2 \ defined by 



-(p) 



^2n T nTn 



l//< 



p=l,2. 



(1) 



For a graphene sheet with a large aspect ratio (width W 
^> length L) we calculate 



Z^^ =0 .85-, r( 2 >=0.87-. (2) 
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Both times are below L/v, as a manifestation of the Hart- 
man effect. 
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Our main interest is in the case of massless electrons, 
but since carriers in graphene may acquire a mass for cer- 
tain substrates [H, , we will include a possible nonzero 
mass term in the calculations. The effect of a mass is only 
important near the Dirac point, so we may set the mass 
to zero in the contact regions, taking the mass profile 
H&{L/2 - |x|) shown in Fig.H 

The electron spin precesses in the y — z plane around 
the magnetic field B = Bx. We assume that the length 
L of the region between the contacts is large compared 
to the length of the contacts themselves, so that we may 
neglect the precession in the contact region and take the 
magnetic field profile B8(L/2 — \x\). 

The Hamiltonian is given by 




-L/2 L/2 



FIG. 2: Profile of the potential V (upper panel), mass fj, 
(middle panel), and magnetic field B (lower panel) along the 
graphene sheet. A nonzero mass is included for the sake of 
generality, but the case [i — is our main interest. The con- 
tact regions | as | > L/2 are modeled by a deep potential well 
(depth Voo 3> hv/L). Spin precession in the contacts is ne- 
glected, so we set B = there. For charge neutrality in the 
region \x\ < L/2 the Fermi energy is lowered to ef — (Dirac 
point). 



H = I a ®H D - \hw B o x ® l ps , (3) 

where I s and I ps are identity matrices in real spin space 
and in pseudospin space, respectively. The Pauli matrix 
cr x in the second term acts on the real spin and accounts 
for the Zeeman energy, with ujb — 9^bB the Larmor 
frequency, (is the Bohr magneton, and g w 2 the gyro- 
magnetic factor. The first term contains the Dirac Hamil- 
tonian, 

H D = V((7 x p x + UyPy) + (T z fl + V, (4) 



II. 



SPIN PRECESSION THROUGH A 
GRAPHENE SHEET 



We study the four-terminal geometry |f4j of Fig. [TJ 
in which spin-up electrons are injected into a graphene 
sheet from ferromagnetic contact f at an elevated voltage 
Vij and drained to ground via three other ferromagnetic 
contacts 2,3,4. The two contacts at the same side of 
the graphene sheet have antiparallel magnetizations. In 
the existing experiments [H, 0] , the contacts at opposite 
sides of the graphene sheet are collinear. This is the 
geometry shown in Fig. [TJ where the magnetizations in 
all four contacts are aligned along the ±y-direction. We 
will consider this case first, and show that it measures 
the time of Eq. {TJ. 

The time is measured if the magnetizations in con- 
tacts 3 and 4 are aligned perpendicularly to those in con- 
tacts 1 and 2 (along the indirection of Fig.[TJ. We defer 
a discussion of that geometry to Sec. IIVI 

Following Ref. Q, the contact regions are modeled by 
a deep potential well, V — —Voo6(\x\ - L/2). (We will 
eventually take the limit Voc — ► oo.) The Fermi energy 
is tuned to the Dirac point sp = in the region between 
the contacts \x\ < L/2. There are therefore no propa- 
gating modes in this region, while the contacts support a 
large number iV^ = h^W/ir of propagating modes (with 
fcoo = Voo/hv the Fermi wave number in the contact re- 
gion). 



for a single valley in graphene (no intervalley scattering). 
The Pauli matrices a x ,u y , er 2 in Hr> act on the pseu- 
dospin (or sublattice) degree of freedom. We neglect the 
coupling between the real spin and the orbit, which is 
weak in graphene. 

We seek the currents /31 and /41 flowing from contact 
f into contacts 3 and 4 separated by a distance L. These 
are determined by the transmittances Tj^ and Tj-f with 
and without spin flip: 



/ 3 i — Tli-GoVi, /41 — T^GqVi 



(5) 



(The conductance quantum Go = 2e 2 /h accounts for a 
two- fold valley degeneracy.) 

For any precessing spin, the probability of a spin-flip 
after a time t is j(ojBt) 2 to second order in B. This 
suggests the definition of an effective transmission time 
t^ 2 \ in terms of the fraction Tj^/(Tjf + T^) = /T^ + 
0(B 4 ) of transmitted electrons that have flipped their 
spin: 



r iT /r TT = i(c B r( 2 )) 2 + o(s 4 ). 



(6) 



Our goal is to calculate this time r 



(2) 
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III. CALCULATION OF THE TRANSMISSION 
TIME FROM SPIN PRECESSION 

The eigenvectors of the Hamiltonian ([3]) corresponding 
to the eigenvalue Sf read 

_ k + iq _ _ k + iq 

Zk ~ EF-V + ti + B' Zk ~ e F -V + f i-B- { ) 

We abbreviate = Hujb/2 and set 7k; to unity (restoring 
units in the final expressions). The wave vectors 

fc* = y/(e F — V± B) 2 — \i 2 — q 2 (10) 
I 



are the longitudinal wave vectors. The wave vector q is 
the transverse wave vector. 

The left spinor in the tensor product in Eqs. ([7]) and © 
represents the state of the real spin and the right spinor 
represents the state of the pseudospin. The superscripts 
+ and — indicate the spin polarization along the x axis: 
the wave functions and 'J' - are eigenstates of a x ® I ps 
with eigenvalues +1 and —1, respectively. 

We solve the scattering problem with potential, mass, 
and magnetic field profiles as shown in Fig. [51 In the 
contact regions |i| > L/2, where V = —Vex, — > oo, we 
have z k + — > 1, z_ k + — > — 1, z k - — > — 1, and z_ k - — > 1. 
We consider a wave incident on the charge-neutral region 
\x\ < L/2 from ferromagnetic contact 1, so with spin up 
along the y-direction. Matching modes at x = ±L/2 wc 
arrive at the following linear equations for reflection and 
transmission amplitudes: 





The amplitudes ru, ri\, £31, and £41 are the reflection 
and transmission amplitudes from contact 1 to contacts 
1,2,3, and 4. Together with the coefficients Ai, A2, A3, 
and A4 we have 8 unknowns, determined by the 8 inde- 
pendent equations contained in Eqs. (TTTj) and (TT2"j) . 
At the Dirac point, that is when ep = 0, we find 



T 
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M 2 = 


4£ 2 k 2 sinh 2 (Lk) 


(13) 


[k 2 - £ 2 + (k 2 + 2? 2 ) cosh (2Lk)] 2 ' 


l<4l| 2 = 


4k 4 cosh 2 (£«;) 


(14) 


[k 2 - £ 2 + (k 2 + B 2 ) cosh (2Lk)] 2 ' 


hl| 2 = 


4B 2 ( K 2 +6 2 )sinh 4 (i/c) 


(15) 




[k 2 - £ 2 + (k 2 + B 2 ) cosh (2Lk)] 2 ' 


kill 2 = 


k 2 (k 2 + £3 2 ) sinh 2 (2Lk) 


(16) 


[k 2 - £ 2 + (k 2 + # 2 ) cosh (2Lk)] 2 ' 



y = W. (Since we assume W 3> £, the choice of bound- 
ary condition does not matter for our results.) The trans- 
verse wave vector is then discretized as q n — 2irn/W, 
where n = 0, ±1, ±2, . . . •|iV 00 numbers the transverse 
modes. The transmittances Tj^ and 7^ j (with and with- 
out spin flip) are defined by the sum over modes of Tj_j 
and T||. For W ^> L and Noo — * 00 the sum over trans- 
mitted modes may be replaced by an integral over q: 



We have abbreviated k = \J q 2 + fi 2 — B 2 . One can verify 
that i? TT +i? iT +T iT +T TT = 1, as it should be. For B = 
(no precession) we recover the transmission and reflection 
probabilities of Refs. @, EH- 

We apply periodic boundary conditions at y = and 



Expanding up to second order in BL = lubL/2v, we 
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(j, [hv/L] 

FIG. 3: Dependence of the transmission time on the mass 
u of the carriers in graphene (for W/L ^> 1). The time is 
below L/v for all fi, as a manifestation of the Hartman effect. 
(While the plot is for r (2) from Eq. ([21]), the time r (1) from 
Eq. (|24p differs only by a few percent.) 




FIG. 4: Dependence of the transmission time for /i = on 
the aspect ratio W/L of the undoped region. Both r' 1 ' and 
r' 2 - 1 are plotted. The limiting values for W/L — > oo are given 
byEq. ©. 



obtain the weak-field transmittances, 



tanh 2 s]u 2 + £ 2 



ttL J (u 2 + £ 2 ) cosh 2 ^/u 2 + £ 2 : 



W 



rr 



du ■ 



J cosh 2 a/u 2 + £ 2 



(17) 



(18) 



/"du . tanh V^+C 2 



7 ^u 2 + £ 2 cosh 2 ^u 2 + £ 



B 2 L 



,2W 



du ■ 



- e2 cosh 2 y/u 2 ->- ' - 
tanh 2 xA^-K 2 



ttL 7 (u 2 + e ) cosh 2 ' 



(20) 



with £ = Lu and m = Lq. 

Comparison with Eq. ^ gives an expression for the 



transmission time 



A 2 ) - 



(L/t 



tanh 2 0/ 2 + £ 5 



OC 

/ 



(/;/ 



(u 2 + e) cosh 2 v/u 2 + £ 2 

1/2 

1 



1/2 



COSh 2 yju 2 + C 2 



(21) 



plotted in Fig. [3l For massless electrons (£ = 0) this 
reduces to 



T ( 2 ) = - 
v 



du 



tanh 2 u 
u 2 cosh 2 u 



1/2 



= 0.87L/v, (22) 



as announced in Eq. @. In the \axge-u limit r^ 2 * 1 — > 
independent of the distance L over which the electrons 
are transmitted. This is the electronic analogue of the 
Hartman effect @, Q . 

These results are for aspect ratios W/L ^> 1, but the 
dependence on the aspect ratio is rather weak, as illus- 
trated in Fig. [4] 



IV. THE CASE OF PERPENDICULARLY 
ALIGNED MAGNETIZATIONS 



We now turn to the case that the magnetization at the 
two ends of the graphene sheet is mutually perpendicu- 
lar, as well as being perpendicular to the magnetic field 
B. Referring to Fig. [TJ we would have the magnetization 
in contacts 1,2 along the ±y-direction and the magne- 
tization in contacts 3, 4 along the indirection (with B 
along x). The transmittances or 7jj defined in Eq. 
([5]) now refer to the transmission of a spin-up in the a y 
basis to a spin-down or spin-up in the a z basis. 

A spin which is initially aligned along the y-direction, 
acquires after a time t a polarization in the z-direction 
given by wst + 0(B 2 ). Analogously to Eq. ©, we now 
define the effective transmission time by 



7j. T _ 7^ 



LO B T 



(1) 



0(B 2 ). 



A very similar calculation as in Sec. IIHI gives 



(23) 



-d) 



(L/v) J du 



du ■ 



tanh \J~v 



^fu 2 + £ 2 cosh 2 ^u 2 + £ 2 
1 



cosh y/u 2 + £ 2 



(24) 



The /x-dependence of r' 1 ' is only a few percent different 
from that of (plotted in Fig. [3]). In the limit £ = 



5 



L/i — > of massless electrons we find 

OO 

m L f , tanhw 7C(3) L 
- y ' = — I au = - — 



L 



-0.85- (25) 
v J u cosh u n z v v 





as announced in Eq. 



Combination of Eqs. (|28|) and ([29)1 shows that 



dq 



tanh p (ivV 



vP{q 2 + fi 2 )P/ 2 cosh 2 (L^q 2 + fj, 2 ) 
-i -l 

dg — — -J— == , (30) 
cosii (LyV + ^i 2 ) 



V. COMPARISON WITH WIGNER-SMITH 
DELAY TIMES 

We wish to derive the relationship (fTJ) between the 
transmission time measured in spin precession and 
the mode-dependent Wigner-Smith delay times r n . By 
definition, the Wigner-Smith delay times are the eigen- 
values of the Wigner-Smith time-delay matrix 



Q 



dS 

dev 



(26) 



constructed from the energy dependent scattering ma- 
trix S. The eigenvalues r„ of Q appear in certain trans- 
port properties [l?], [3 , but they are usually not directly 
measurable. For example, the thermopower of a single- 
channel conductor depends on the difference n — T2 of the 
two eigenvalues of Q, as well as on the eigenvectors. It 
is therefore not obvious a priori that can be related 
to the t„'s. 

Since we seek the delay times in the limit of zero mag- 
netic field, we can consider a simpler scattering problem 
than in the previous section, namely transmission of spin- 
less electrons through a graphene sheet with the mass 
and potential profile shown in Fig. [2j In this case the 
scattering matrix is given by [7fl 



I 



k cos kL — iep sin kL 

\-q-i^)smkL k . 

k (q — ifi) sin kL J ' ^ ' 



where k — y/ej? — A* 2 — q 2 . 

The general energy-dependent expression for Q is 
lengthy, but at the Dirac point it simplifies to 



Q = t®{0 lj< T(g) = vV q 2 + , 2 ' (28) 

So for each mode n there is a single doubly degenerate 
Wigner-Smith delay time t„ = r(q — 2nn/W). The 
mode-dependent transmission probability at the Dirac 
point is T n = T(q = 2-ku/W) with 



T(q) = 



1 



cosh 2 (L^/q 2 + /i 2 ) 



(29) 



where we have replaced the sum over modes by an in- 
tegration over wave vectors (appropriate for W/L ^> 1). 
Comparison with the expressions (|2ip and (|24|) for 
and proves the identity (TTJ) of the transmission time 
measured in spin precession and the weighted average of 
the mode-dependent Wigner-Smith delay times. 



VI. CONCLUSION 

In conclusion, we have shown how spin precession in 
graphene may reveal an unusual dynamical aspect of bal- 
listic quantum transport at the Dirac point. In a clean 
charge-neutral graphene sheet of length L, the transmis- 
sion is via evanescent rather than propagating waves. 
While for propagating waves the transmission time is 
bounded by r > L/v, evanescent waves have no well- 
defined velocity and can show a shorter r in a precession 
measurement. This is the electronic analogue of the Hart- 
man effect from optics [H, [l(| • Our result ^ for massless 
electrons is not much below r = L/v, but it does provide 
an unambiguous demonstration of this apparent superlu- 
minality. 

From a conceptual point of view, our analysis demon- 
strates, firstly, that the pseudodiffusive aspects of ballis- 
tic transmission at the Dirac point (as observed in con- 
ductance and shot noise [H, 0j H3), are restricted to 
static properties. The dynamics is not diffusive in any 
sense (no L 2 scaling of r). Secondly, our analysis demon- 
strates via the relation (TTJ) that the Wigner-Smith delay 
times are directly observable through spin precession at 
the Dirac point. 

We finally notice a qualitative difference between spin 
precession in a tunnel barrier and spin precession at the 
Dirac point. As pointed out by Biittiker [22j], the spin of 
a tunneling electron not only precesses in the y — z plane 
perpendicular to B, but in addition aligns itself along 
the magnetic field. The rotation of the spin out of the 
y — z plane (dominant in a tunnel barrier, but ignored in 
the Larmor clock @, Q) appears because of a difference 
in tunnel probabilities for spins parallel or antiparallel to 
B. No such out-of-plane rotation appears at the Dirac 
point, due to the fact that the energy-dependent trans- 
mission probabilities are extremal at zero energy. The 
spin precession geometry analyzed in this work is there- 
fore particularly close to the original concept of a Larmor 
clock. 
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